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1. Introduction and Motivation 

The investigation of the gluon propagator in the Landau gauge, with particular emphasis on 
the deep infrared region, has been a hot topic in the recent years, with many contributions coming 
from lattice QCD simulations and the solutions of the Dyson-Schwinger equations. 

The QCD simulations for studying the gluon propagator in the Landau gauge with the SU(3) 
and SU(2) gauge groups, typically, use the Wilson action, which approaches the continuous ac- 
tion up to corrections of order €?(a 2 ). Furthermore, the gauge procedure ends up satisfying the 
lattice version of the Landau gauge condition d ■ A = 0. Again, the lattice gauge fixing condition 
reproduces the continuum condition up to corrections of order ff(a 2 ). 

In order to access the deep infrared properties of the gluon propagator, simulations have been 
performed on large volumes with a lattice spacing a ~ 0.2 fm or, equivalently, aT x ~ 1 GeV which, 
for SU(3) simulations, corresponds to a j8 =5.7. From the point of view of the lattice simulations, 
such a small j8 stands at the beginning of the perturbative scaling window and, therefore, one should 
expect a relatively large contribution from the finite lattice spacing effects. 

In §\§ the finite lattice size and anisotropic effects were studied using a simulation with j8 = 6.0, 
i.e. a = 0.10 fm. Methods to minimize these artifacts were also investigated. These methods, the 
cylindrical and conic cuts, are now widely used. This study was repeated in [Q] together with 
the investigation of finite volume effects by performing simulations for j8 = 6.0 and j8 =6.2, i.e. 
a = 0.075 fm. If the finite lattice spacing effects seem to be under control, the finite volume ef- 
fects show up mainly in the infrared region, with the zero momentum gluon propagator decreasing 
with the lattice volume. Similar conclusions were obtained using larger lattices volumes - see, for 
example, [||, [|, 26, |6fl. If this effect is clearly seen in many lattice simulations, the extrapolation 



to the infinite volume is not so unanimous. In particular, depending on how you extrapolate to 
the infinite volume [f7|], one gets D(0) = 0, as predicted in [|8[ |9|], or D(0) ^ 0. Generalizations 
of the Zwanziger work do not require necessarily a vanishing zero momentum gluon propagator 
[10, 11, |12|]. Lattice simulations at finite volume, modulus the extrapolation to infinite volume, are 



in good agreement with such scenarios [13]. On the other hand, looking at quarks models derived 



from QCD that incorporate the propagation of the gluon [14], it seems that phenomenology does 
not care about the value of D(0). Given that the zero momentum gluon propagator is connected 
with gluon confinement, it is a fundamental quantity and, in this sense, it would be nice to be able 
to have an estimate for D(0) from lattice simulations free of lattice artifacts. 

In section ||, we review the large volume simulations, i.e. having at least V = (5.8 fm) 4 , for the 
gluon propagator performed by the authors and discuss the possible values for D(0). 

Besides the infrared behavior of the gluon propagator, we would like also to discuss the mass 
scales in the gluon propagator. The lagrangean of QCD has no mass scale. However, dimensional 
transmutation introduces a scale, Aqcd, in association with the ultraviolet properties of the theory. 
Aqcd is used to separate phenomena which belong to the non-perturbative regime of the theory, 
i.e. where the relevant momenta are small compared to Aqcd, or to the perturbative regime. 

From the point of view of the gluon propagator, if one wants to parameterize its functional form 



more mass scales are required, being it as an effective gluon mass [15] or masses [16], the Gribov 



mass parameter m, W or via the introduction of masses which are related with condensates required 



to define a non-perturbative Green's function generator for QCD [10, 11, 12]. For example, in [ 13] 
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Figure 1 : Renormalized gluon propagator. The numbers in the picture are the values for D (0) for the /3 = 5 .7 
simulation using a 44 4 lattice, j8 = 6.0 with a 80 4 lattice and J3 = 6.2 for a 80 4 lattice. 



the non-perturbative part of the lattice gluon propagator was described by 

D(q 2 ) 



q 2 +M 2 



q 4 +Mjq 2 +M 4 



(1.1) 



with the mass scales being M\ = 2.15(13) GeV 2 , M\ = 0.337(47) GeV 2 and M 2 = 0.51 GeV 2 . 
The precise meaning of these mass scales depends on the theoretical scenario under discussion. 
The gluon propagator as given by equation ( |1.1| ) was used in [17] to compute the lowest mass 



glueball states within the so-called Refined Gribov-Zwanziger theory. Furthermore, in [18, 19] a 
Gribov-type propagator 



D(q 2 ) 



(1.2) 



q 4 +M 4 

was used to investigated the QCD operators which can correspond to relevant physical degrees of 
freedom, i.e. those operators whose propagator have a Kallen-Lehmann representation. The Gribov 
type of propagator was also explored in [|(]] to develop a theoretical scenario to explain gluon 
confinement. Further, non-perturbative solutions of gluon-ghost Dyson-Schwinger equations also 



allow for a dynamical generated gluon mass [ |21| , [22| , |23[ ] which, again, require other mass scales 
than Aqcd- 

In section |3[ we try to define the minimum set of mass scales which can describe the lattice 
gluon propagator over the full range of momenta using two possible scenarios. For the two sit- 
uations considered, we find that the propagator requires the introduction of four mass scales. It 
remains to see what are the relevance of these masses, if any, to physical phenomena. 

2. Gluon Propagator from Large Volume Simulations 

At Coimbra we have been computing the Landau gauge gluon propagator for several lattice 
spacings and volumes. Our goal is to perform a combined analysis of the finite lattice spacing 
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Figure 2: Renormalized gluon dressing function. 



and volume effects. We have just beginning such an investigation and the results will be reported 
elsewhere. Here, we will focus on our large volume simulations. 

Figure [j] shows the renormalized gluon propagator for our largest lattice volumes for several 
/3 values. The corresponding values for D(0) are also reported. Definitions and the details on the 
computation of the propagator can be found in [j24p. 

The propagators define a unique curve for momenta q > 0.7 GeV. For smaller momenta, the 
propagator data for the finest lattice, which has the smallest physical volume, is below the corre- 
sponding data points for the j3 = 5.7 and j8 = 6.0 simulations. This is certainly a finite volume 
effect which deserves to be further investigated. A closer look at the so called gluon dressing 
function q 2 D(q 2 ), see figure ||, shows small but finite volume effects. 

For smaller lattice volumes, the gluon propagator follows the behavior observed in figure []]. 
In all the simulations, no turnover on the gluon propagator at small momenta is observed. In what 
concerns the value for D(0), the values reported are, within two standard deviations, compatible for 
all simulations. A combined value, averaging taking into account the statistical errors, gives D(0) = 
8.536(70) GeV" 2 . This value is in good agreement with the estimate of [ |13| ] where D(0) = 8.3(5) 
GeV 2 and slightly above the result of the large volume simulations performed by the Berlin- 
Moscow-Adelaide group p5|] where D(0) = 7.53(19) GeV" 2 . 



Can these results be an indication that D(0) is finite and non- vanishing in the infinite volume 
limit? To our mind, that depends on the way you want to look the gluon data. If one takes the 
propagator at its face value, then the lattice simulations points clearly towards a finite and non- 
vanishing D(0). However, one should not forget that in the deep infrared region, one is relying on 
on-axis momenta and one is banging on the lattice walls. Anyway, from the simulations performed 
so far, the results look rather stable. However, if one looks at some details of the calculation, there 
are some holes which should be investigated, in particular, the question of zero modes which can 
give a large contribution to the deep infrared. 
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Figure 3: Ratios of the gluon propagator. See [pjj for details. In the left hand side, the plots shows the 
ratios as directly obtained for the gluon propagator. In the right side, the ratios were shifted in order to agree 
for the highest R q [n] . 



In [26] we have looked at the compatibility of the lattice gluon and ghost propagators with 
a pure power law behavior in the deep infrared region using large asymmetric lattices. We will 
not reproduce the work here but will resume its main points. In the above cited work, we have 
investigated ratios of propagators as a function of ratios of momenta. The ratios of propagators are 
renormalization independent and, if the propagators are described by a pure power law, are univer- 
sal functions. Furthermore, the advantage of the ratios being that they suppress finite volume/lattice 



spacing corrections. The results reported in [ |26| ] show that while the gluon propagator seems to 
follow a pure power law in the deep infrared, the ghost propagator doesn't. Moreover, at least for 
the gluon propagator, the results suggest a way to correct the propagator data. If this is taken into 
account as described in the appendix of [pq], not only the corrected propagators become closer to 
each other, but also the corrected propagator decreases for momenta below 400 MeV. 

Figure || shows the ratios for the gluon propagator computed for several lattice simulations. In 
the left hand side, the plot shows the ratios as measured from the gluon propagator. The differences 
between the different data can be understood as due to finite volume effects. In the right hand side, 
the propagator data were shifted to match the point at the highest R q [n] , which corresponds to the 
smallest non- vanishing momenta available for each lattice. From the right part of figure [| one can 
claim that the lattice propagator defines a universal curve and, in this sense, claim that the gluon 
propagator is compatible with a pure power law in the deep infrared. According to this analysis, 
the deep infrared propagator is described by 

2\2k 



with k = 0.528(14) 



(2.1) 



and D(0) vanishes. Further details can be found in [|27|. 



From the point of view of the phenomenology, the results of [|14|], although obtained within a 
quark model, suggest that it would be quite hard to get some insight from phenomenology on the 
"correct" value for D(0) and distinguish, in this way, what value one should consider. 



3. Mass Scales from the Lattice Gluon Propagator 

We start our discussion of the mass scales for the lattice gluon propagator assuming that it can 
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Figure 4: The infrared gluon propagator as a massive type propagator. 



be described by a massive type propagator [|I5|], i.e. that 

Z 



D{q 2 ) 



q 2 +M 2 



(3.1) 



The lattice data is fitted, in the range q £ [0,g OT(rc ], to equation (|37TJ). Figure || resumes the outcome 
of fits. 

From fitting the gluon propagator computed with the largest physical volume, i.e. the simula- 
tion using j8 = 6.0 and an 80 4 lattice, it follows that the lattice data can be described by equation 
( |3~T| ) up to momenta q = 482 MeV, which has the smallest % 2 /d.o.f, with M = 641(10) MeV and 
Z = 3.99(13). The corresponding % 2 jd.o.f . being 1.15. This parameters give aD(0) = 10.58(41) 
GeV~ 2 , which is slightly above the values claimed in the previous section. 

The full range of momenta cannot be described by a massive type propagator, unless we al- 
low Z and M to become functions of the momenta. The computation of Z(q 2 ) and M 2 {q 2 ) were 
investigated in [15]. Performing a slide window analysis, it was found that Z{q 2 ) can be described 
by 



Z(q Z 



[A + ln (q 2 +m 2 )] : 



(3.2) 



0.4 and m\ ~ 1.57 



where 7 is the anomalous gluon dimension, and the fits give Zq ~ 1 .0, A • 
GeV 2 for a % 2 jd.o.f. = 1.8. The mass squared is well described for momenta above 1 GeV by 



') = CQ + ciq 2 Inq 2 



(3.3) 



with q given in GeV, cq ~ —0,3 GeV 2 , c\ = —0.2 for a x 2 jd.o.f . = 1.8 and, for momenta below 1 
GeV, is given by 

M 2 (q 2 )= 0.58(4)- 1.24(4) q 2 ; (3.4) 
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Figure 5: The lattice gluon dressing function q D(q ) as a sum of multi massive type propagators. The 
various curves correspond to various N. The legend includes the mass scales corresponding to each fit. 



the x 2 /d.o.f '. = 1.8. Note that (3.4) is the expression of [21] for the dynamical generated gluon 
mass, where a non-perturbative solution of the gluon-ghost Dyson-Schwinger equations was inves- 
tigated. 

Another possibility to describe the gluon propagator is to assume that D(q 2 ) is given by a sum 
of massive type propagators, i.e. 



N 



D(q 2 



y _±! 

u q 2 +M 2 



(3-5) 



i=\ 



Note that equation ( p.5| ) is a simplified Kallen-Lehmann representation for the propagator. The 
figure shows that by increasing N one is able to described the lattice data over a wider range of 
momenta. The various fits give 



0.57(^{3.535(64),0.5907(86)}J 1.4 
1.52({17(3),0.797(17)}{-17(3), 1.035(31)}) 1.5 
6.46({31(6),0.851(16)}{-52(11),1.062(26)}{22(9),i.257(40)}J 1.6 
7.77 f{33(9),0.900(26)}{-54(12),1.163(49)}{33(14),1.65(12)}{-ll(ll), 2. 11(24)}) 1.1 



where the results were written as q max ^{Z,-,M,},/ = 1 • • -Nj % 2 /d.o.f. , with q and M expressed 
in GeV. The last line corresponds to the full range of momenta for the largest physical volume that 
we have simulated. We would like to call the reader attention to the stability of the numbers, when 
increasing N and the fitting range. Furthermore, note that the lattice gluon propagator requires 



some negative Z,. Given that (3.5) can be viewed as a Kallen-Lehmann representation, negative Z 
means that the gluon propagator is not an asymptotic state of the 5-matrix. 

The above investigation shows that the lattice gluon propagator requires four mass scales. Note 
that the values for the masses are within typical hadronic masses. 
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Figure 6: The lattice gluon dressing function q~D(q~) as a sum of multi massive Gribov type propagators. 
The various curves correspond to various N. The legend includes the mass scales corresponding to each fit. 



As a final comment, we would like to address whether the gluon propagator can be described 
by a multiple mass Gribov type formula, i.e. 



N 



(3.6) 



A similar investigation as for expression ( |3.5| ) shows that when N = 1 the Gribov formula does not 
match the lattice data, i.e. for /V = 1 it turns out that % 2 jd.o.f . is always above 2. However, if 
one considers higher N, then the multiple mass Gribov type propagator ( |3T6| ) is compatible with 
the lattice data in an appropriate range of momenta. Following the same notation as before, the fits 
give 

0.90( {0.430(19), 0.1639(38)}, {1.692(17), 0.4343(31)}) 1.6 



3.60 (^{0.487(19), 0.1724(36)}, {1.766(17), 0.4588(37)}, {-1.3412(83), 1.505(15) M 1.6 
7.77 ( {0.504(20), 0.1748(37)}, {1.773(17), 0.4647(43)}, 



{-1.306(12), 1.450(12)}, {-0.251(12), 4.48(19) }J 1.3 

and can be seen in figure |6| together with the lattice data. Again, the description of the lattice data 
over the full momentum range requires four mass scales. 

In what concerns the value for 0(0), the multiple massive type propagator gives a 0(0) ~ 10 
GeV 2 (with large statistical errors), while the multiple Gribov type propagator requires 0(0) = 0. 
The Gribov-type propagator also provides an estimate where the turn over can be observed. Indeed, 
the first derivative of D(q 2 ) vanishes for q = 190 MeV, with the derivative being negative for 
momenta above 0.19 GeV and positive below. The renormalized gluon propagator together with 
the fit to (3.6) can be seen in figure [7[ 
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Figure 7: The lattice gluon propagator as a sum of multi massive Gribov type propagators. 

4. Final Comments 

In this article we review the calculation of the Landau gauge gluon propagator by 00 and 
PJS performed in Coimbra using lattice QCD simulations. Although, the lattice data seems to 
point towards a finite and non-vanishing D(0) around 8 GeV -2 , there are still some open questions 
which deserve to be investigated. In particular, why the ratios of propagators point towards a power 
law behavior in the infrared region and a D(0) =0 needs to be understood. 

In the second part of the article, the modeling of the lattice data using massive type propaga- 
tors and Gribov type propagators is discussed. It is shown that the full range of momenta can be 
described using four different mass scales. On the other side, assuming a mass dependent gluon 
mass, its functional dependence is sketched. In particular, it is shown that the lattice data is compat- 
ible with the parameterization used in the non-perturbative investigation of the Dyson-Schwinger 
equations. In what concerns the infrared gluon propagator, we show that up to momenta around 
500 MeV, the gluon propagator is well described by a massive type propagator, with an effective 
gluon mass around 640 MeV. 
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